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5_i ■ We discuss how concepts such as geodesic length and the volume of space-time can appear 



in 2d topological gravity. We then construct a detailed mapping between the reduced Her- 
mitian matrix model and 2d topological gravity at genus zero. This leads to a complete 
solution of the counting problem for planar graphs with vertices of even coordination num- 
ber. The connection between multi-critical matrix models and multi-critical topological 
gravity at genus zero is studied in some detail. 
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1 Introduction 

Despite substantial progress in our understanding of two-dimensional quantum gravity we 
are still left with two separate theories: two-dimensional topological gravity and Liouville 
theory. Topological gravity is a theory where concepts such as geodesic distance make no 
sense. It can be viewed as a theory on the moduli space of punctured Riemann surfaces, 
and it is defined in such a way that it allows us to calculate intersection indices related to 
the moduli spaces. In Liouville theory geodesic distance plays a key role: Hartle-Hawking 
wave functions are defined in terms of the geodesic length of the boundaries |I|, |2[ and it is 
possible to consider the quantum theory of universes of fixed space-time volume, as well as 
correlation functions depending on the geodesic distance between two punctures Q (see 
also references j|, |j for a discussion of more general correlation functions). Nevertheless, 
the two theories seem to be identical. The various proofs of the equivalence have so far 
been rather indirect || [/J. Typically, generalized Schwinger-Dyson equations have been 
derived for the two theories and have been shown to be identical by suitable identifications 
of variables. By such identifications the concepts of length and area enter into topological 
gravity without any obvious interpretation. As an example let us mention the following 
situation: in topological gravity we can calculate the expectation values (a n ) of certain 
generalized puncture operators and we can define the generating function for this set of 
observables by 

pn 

(w(i)) = £ - K> • (i) 

n n - 

In topological gravity £ has no interpretation, beyond that of an indeterminate which 
defines the generating function. However, {W{£)) agrees with the corresponding Hartle- 
Hawking wave functional calculated by the use of dynamical triangulations. In Liouville 
theory £ has from first principles an interpretation as the length of a macroscopic boundary 
of the two-dimensional universe. The same arguments apply to the total space-time 
volume A, which is defined from first principles if one uses Liouville theory, while it has 
no natural definition within the framework of topological gravity. 

The rest of this paper is organized as follows: in section y we discuss how concepts such 
as the volume of space-time and the length of boundaries can be introduced in topological 
gravity. Since these concepts are natural in topological quantum gravity it should be 
possible to make a detailed mapping to Liouville theory. This is most conveniently done 
if we use the dynamically triangulated version of Liouville theory. In order to define 
the mapping we introduce in section |3| a contour integral representation of the partition 
function of topological gravity. A related integral representation for the Hermitian matrix 
model allows us in section ^| to perform a simple and detailed identification between the 
correlation functions in topological gravity and those in quantum gravity defined via 
dynamical triangulations. As a side result of this identification we use topological gravity 



results to solve completely the counting problem for planar graphs, whose vertices all 
have even coordination numbers. In section |5| we apply the mapping between topological 
gravity and matrix models to the case of multi-critical models. Section ^ contains our 
conclusions. 

2 Lagrangian Aspects of Topological Gravity 

Various starting points are possible in the case of topological gravity. One approach is 
to start directly from geometric considerations involving the intersection indices on punc- 
tured Riemann surfaces H 9j. This approach was completed and made mathematically 



rigorous by the seminal work of Kontsevich [1C]. From the remarkable properties of the 



Kontsevich integral it follows that the partition function of topological gravity is a tau- 
function for the Korteweg-de Vries equation. More precisely, the partition function is 
given by the following matrix integral 

where the integration is over N x N Hermitian matrices <j). The matrix M is assumed to 
be symmetric and positive definite, and one can show that Z^{M) only depends on 

t k = - {2k ~ 1] " TrM-( 2fc+1 ), k>0, (3) 

where by convention (—1)!! = 1. For N — ► oo the u- function of the KdV-equation is then 

u({t t }) = -^logZ({ tl }). (4) 



'o 



One can write (for N — > oo): 



log z({ ti }) = Upj2 *i°i ) = E E • • • K vr • • •> n -fr . ( 5 ) 

\ j=0 I n =0ni=0 j=0 \ ■?'/ 

This expression has a large N expansion which is at the same time a genus expansion and it 
allows us in principle to calculate the numbers (ctq ^ 1 ■•■)„, which have the interpretation 
as intersection indices for each genus. 

A number of proofs || exist of the equivalence between topological gravity formulated 
via (Q) and Liouville gravity formulated via the use of dynamical triangulations. As al- 
ready mentioned the proofs of the equivalence between the theory defined by (g) and the 
theory of gravity defined via the matrix models corresponding to dynamical triangula- 
tions^] are rather indirect and leave no hints of how to define observables associated with 

4 We view here Liouville gravity and simplicial gravity as representing the same kind of quantum gravity, 
since observables referring directly to geometric objects (Hartle-Hawking wave functionals etc.) can be 
defined in both theories, and since all calculations in the two theories agree so far. 



geodesic distances in topological gravity. For this purpose it is more convenient to turn to 
the original formulation [11| of topological gravity, which originated from Witten's general 



approach to topological field theories [12]. Here the starting point for a topological theory 



is 



C = 0. (6) 



The non-trivial content of (|6J) comes by defining the space of variables and performing a 
corresponding gauge fixing of the local symmetries. In the case of gravity one can choose 
as the field variable the metric g^ u . The Lagrangian ([]) is invariant under the most 
general variation 

Sg^v = Kfjtu- (7) 

It is common to factor out explicitly the invariance under diffeomorphisms. If £^ denotes 
an infinitesimal diffeomorphism and r„ a general 2x2 matrix, one can write 

$9vw = £fi-,v + &/;/* ~ (r^g\u + r^gxfj), (8) 

or, relating more closely to non-critical string theory, 

&911V = £/x;i> + &;/x ~ ^9/i^X + 9^, (9) 



in decreasing order of generality (see |y|, [Tj] for a discussion of the possible representations 



of (pi)). Gauge fixing according to either (pi) or (pi) amounts to defining (||) to mean 

£ GF = s ■ *, (10) 

where s is the BRST operator for the gauge fixed symmetries. For both gauge fixings (JSj) 
and (0) the split into diffeomorphisms and a residual shift is not unique, since there exist 
transformations different from the identity which leave g^ u invariant. This leads to the 
well known ghosts for ghosts (super-ghosts) in C GF . In the following we will use (|8]) where 
the corresponding field content of the gauge fixed Lagrangian ((TTJ) is the two-dimensional 
metric g^ u , the gravitino field ip^ (the ghost field of the shift), the ghost field c M of 
diffeomorphisms and the super-ghost field 7^. In this representation, the action of the 
BRST operator s is as follows W\ 



S ' 9^iu — ■£ c g^iu 1 Ipfiu > S • C — 2 ■*- c C + T > /-.-.n 

S ■ W\iv * c W\iv *"y 9\nv 1 8 ' 7 £ c 7 > 

where £ c and £ 7 denote the action of infinitesimal diffeomorphisms with parameters &" 
and 7^, and where ip^ u is related to the ghost wj 1 for the matrix rj 1 by 

Vv = ~Vsv ~ wJg-yiJ.- ( 12 ) 



One obtains the observables starting from the Euler 2-form [15] 



cr( 2 ) = — y/gR e^dx? A dx u , (13) 

which gives rise to a 1-form and a 0-form via the descent equations: 

s ■ o^ = d a^\ s- a« =da<®, s • ^ = 0. (14) 

In this equation s is the BRST operator and d the exterior derivative. Explicitly one gets 



0.(0) _ V9£y.is 

(1) \J~9 E iiv 

4tt 



1 -c^r + c »D p (r p - g vp r„) + d»y - \r P r p ] , (15) 

(?R + D p {r p -g vp r„)\dx^ (16) 



where D p denotes the covariant derivative. One can build an infinite set of observables 
from a' ' since the super-ghost 7^ is commutative: 

4 o) - M o) ) n (17) 



The expectation values of observables in topological gravity are now defined by func- 
tional integrals like 

/nswWwe-^n^^), (is) 

where $ denotes a suitable collection of fields (see below) and the Lagrangian C associated 
with the action S[$] is BRST exact. The zero-forms a nk sit at the points P^ on the 
world-sheet manifold and the super-gauge transformations should leave these points fixed. 
Correspondingly, the fields c and 7 must vanish at these points and this implies that we 
consider punctured Riemann surfaces. In this way the Riemann surfaces considered in 
the path integral become a function of the observables. This is in accordance with the 
interpretation of the "observables" Tr <fi x • • • TV <ft k in the Hermitian matrix model. 
The matrix integral corresponding to these observables is equivalent to a summation over 
(generalized) triangulations with k boundaries of length 2l\, . . . , 21^, i.e. for each matrix 
observable we consider different triangulated surfaces. For fixed l\,. . . ,1/- the continuum 
lengths of the boundaries go to zero and the loops become punctures. In later sections 
we will make the correspondence between the cr nk 's and the matrix model observables 
precise. 



In order to perform the functional integral (18) one usually chooses a background 
metric g° and background gravitino field i/)°, which depend on the moduli of the punctured 
Riemann surfaces, and after successive integration over moduli all explicit reference to 
the points P^ disappears [O]. In the following we will explicitly suppress all dependence 



on moduli as they play no role in our arguments. The background metric and background 
gravitino are introduced in a BRST exact way via Lagrange multiplier anti-ghosts b^ v and 
Pftn, and their corresponding Nakanishi-Lautrup fields d^ v and 5^ v . The BRST algebra 
extends trivially to these new fields 

s .f3^ = 5^, a- 6^ = 0. (19) 

The gauge fixed action is 

S GF = / d 2 XyfgC G¥ = d 2 x^/gs-^/ 

d 2 x^s- (b^(g^ - g%) + /3^(W " <j) • (20) 



It is a remarkable fact that the functional integral (18) allows the extraction of almost 
complete information about the intersection indices on the moduli space of punctured 
Riemann surfaces. This is done by deriving Ward identities between various correlators 
and by showing that these Ward identities allow a recursive determination of expectation 
values of the a„'s [16], O ]. 

Contrary to the Kontsevich integral formulation this definition of topological gravity 
operates with space-time, but offers a priori no concept of space-time distance. Some- 
times it is even said that there should be no such concept in quantum gravity, since by 
reparameterization invariance, we only know if two points are separated or coincide and 
that this is the origin of the so-called contact terms. This is in contrast to the approach 
to quantum gravity via Liouville theory, or, at the constructive level, via dynamical tri- 
angulations. To be more explicit we can in Liouville theory define the partition function 
for a fixed space-time volume by 

Z W (A) =Jv[g] 5 ^Jd 2 ^-A^j , (21) 

and a Hartle-Hawking wave functional corresponding to a fixed boundary length £ and 
fixed space-time volume A by: 

W hV (£,A)=Jv[g]5^Jd 2 ^-^) sUds-iY (22) 

Similarly we can define the partition function of quantum universes with two marked 
points separated by a geodesic distance R and with space-time volume A by 

G LV (R,A) = Jv[g\ 5 ^Jd 2 ^-A^j J ' d 2 '^ J ' d 2 £ ' J~g 5 (D 9 (£, £') - R) , (23) 

where -D 9 (£,£') is the geodesic distance between £ and £' with respect to the metric g. In 
these equations the functional integration is over all equivalence classes of metrics. In a 

6 



certain sense Z hV (A), W hV (£,A) and G LV (R, A) are "more topological" than topological 
gravity itself in the spirit of being a functional integral with C = 0, since no action is 
ever used in (^l|)-(p3D. However, at the same time they have explicit reference to metrical 
properties of the quantum space-time. In fact G LV (R, A) is a perfect probe of the fractal 
structure of space-time [|| fj. 

We can also construct a partition function for topological gravity which corresponds 
to a fixed space-time volume. To do this we introduce the appropriate delta functions as 
a global, BRST exact term in the action, 

Z[A ] = J V^dxdfi exp | -S GF -s- X (J sfg ~ 4,) } , (24) 

by use of Lagrange multipliers forming a constant anti-ghost multiplet (x, A*), which trans- 
forms like 

s ■ x = M> s- n = Q. (25) 



In contrast to equation (|2l]) this also includes an additional delta function in the gravitino 
field, induced by super symmetry. We must demand that the area of the surface is compat- 
ible with the chosen background metric, that is, Aq = J* s y/g^. By Laplace transformation 



of (24) one finds a new topological action with a BRST exact "cosmological term", 

Z[n] = J V^d x exp |-5 GF - s ■ (x J^ y/gj } • (26) 

If the surface £ has boundaries it is also possible to introduce explicitly a "boundary 
cosmological term" living on the boundary of S fl8fl , and by an inverse Laplace transfor- 
mation like (^4|) one can fix the length of the boundary and in this way obtain a topological 
gravity equivalent of (^2|)p|. 

The above considerations show that one can indeed define topological gravity on a 

restricted moduli space of metrics where for instance the area A has a fixed value. In 

this way there is no conflict between the appearance of observables that refer to global 

metrical properties, such as the total area or the total length of a boundary, and the 

concept of topological gravity. We now proceed to establish a precise mapping between 

the operators a n as they appear in topological gravity and the corresponding operators 

in the language of dynamical triangulations or matrix models. 

5 Strictly speaking it has not been possible to calculate W LV (£, A) and G LV (R, A) entirely within the 
continuum framework of Liouville theory. However, they can be calculated using dynamical triangulations 

ail- 

6 We believe that it is possible to define a topological gravity partition function corresponding to (|23|), 
too, but since we do not know how to calculate the corresponding partition function within the framework 
of topological gravity as a function of geodesic distance, we have not pursued this question further. 



3 Topological Gravity at Genus Zero 

For the purpose of establishing this required map between topological gravity and matrix 
models we will consider genus zero surfaces only. Further it is simplest to start from the 
formulation of topological gravity found in M . 

3.1 Contour Integral Representation 

In reference ||] Witten considers 2d topological gravity with a perturbed action. The 
theory contains a set of basic operators, {cr n }, where n is an integer with n > 0, and 
the perturbation is parameterized by a set of coupling constants, {t n }. The unperturbed 
correlation functions will be written with the subscript V = 0; below we introduce a 
potential V(z), which depends on the perturbed action, and ^ = corresponds to the 
unperturbed case. Then the partition function, (1)„, for the perturbed theory at genus g 
is given in terms of the unperturbed correlation functions for the same genus, 



(l)v = ( exp 



E l 3 a 3 
3=0 



(27) 



v=o 

OO OO / OO J. n 3 \ 

= E E- n:M<M n >i) ni (^r---> F=0 . (28) 

no=0ni=0 \j=o' L rJ 

For genus zero the unperturbed correlation functions can be easily calculated from the 
puncture and dilaton equations |jl| . With the normalization of a n used in reference [|| , 
we have, 

(n°*) =(e*) 1 *f s +E(*- 1 )). ( 2g ) 

where S(n) = S n Q. The correlation function of a product of observables for unperturbed 
topological gravity vanishes unless the total ghost number of the observables equals the 
dimension of the moduli space; the delta function reflects this fact. At higher genus the 
unperturbed correlation functions are much more complicated. The rest of this paper is 
almost exclusively concerned with genus zero surfaces and all the correlation functions 
given are for genus zero, unless it is stated otherwise. 

Now we will define a potential V(z) for the theory given in equation (f27|), 

OO 

V(z) = Y,t3Z J , (3°) 

3=0 

where we have introduced a new complex variable z. We will write to as e and the 
puncture operator ctq as P, when it is convenient to do so. Then we can rewrite the 



partition function for the perturbed theory at genus zero as a contour integral. We have 



(1), 



E E 

no=0 ni=0 

d_ 
d\ 



n 



t; 



rij\ 



J2 kn kV 6[3 + ^T(k-l)n k 



-1 oo oo oo /+"J 

e E-n m^- 1 ^ 

n =0ni=0 j=0 \ 1' 



A=0 
cocff. of z~ " 



where 



1 



d_ 



-1 oo 



E 

n=0 



d_ 



(31) 
(32) 

(33) 



The derivatives with respect to A, give the factor of (J2 fcn&)! an d keeping only the coef- 
ficient of z~ 3 enforces the delta function. Thus, 



(1), 



1 



£H£M>: 



(34) 



J A=0 



where we are integrating around some suitable contour encircling the origin. Changing 
variables z — ► z/A gives 

dAj / 2tt« A 3 V : , , 

-I A — U 

Expanding the exponential in powers of A and performing the differentiation with respect 
to A gives 



(1), 



1 



~ L f dz z 2 (\ , 
exp -V{z, 



(35) 



(1), 



E 

r=3 



dz 



(Viz) 



1 



2tti \ z J r(r — l)(r — 2) 
It is convenient to rescale V by a factor of (3; then after V(z) — ► (3V{z) we have 



^3_ 

d/3 3 



(1) 



,«- 



dz /FM 1 
2^ I T" J 1-pV/z' 



Integrating with respect to /3 three times and then setting (3 = 1 gives 



(1), 



c/z 
2~7ri 



(z - y) 2 log 



z-F 



(36) 



(37) 



(38) 



a 



This formula, which gives the genus zero partition function for 2d topological gravity with 
an arbitrary polynomial perturbation V(z), will be used extensively throughout the rest 
of the paper. For a suitable choice of the parameters {t n }, the integrand has a cut on 
the positive real axis from z = to some point, which we will call z = u. This point is 
defined by the equation u = V(u), that is, 



+ E^ ? 

fc=i 



(39) 



which is just the well-known string equation for genus zero ||. The contour C encircles 
this cut, across which the integrand has a discontinuity of —2tti. It will be convenient to 
assume that e > 0. 



3.2 Multi-Critical Models 

In this section we will calculate various correlation functions for the multi-critical models. 
The simplest choice of potential for the k-th multi-critical model is V = e + z — z , that 
is to 7^ 0) t\ = 1 an d tk = — 1. Note however that this choice of potential is not unique 
and in section |B| we will show that the usual Hermitian one-matrix multi-critical models 
correspond to a different choice of V{z). The scaling behaviour for the k-th multi-critical 
model is however independent of the particular choice of V(z). 

From (p8|) we have that 



c v 



,k 



W=-2 <?— \ Z " e ) lo g(^)' ( 4 °) 



for the /c-th multi-critical model. There is a cut on the real axis in the interval 
which is encircled by the contour C . Thus, 



0,e l / k 



(1) = - / dx (x k - e) = \ ,-, (41) 



2 Jo V ) -(2 + I)(l + I) 

which agrees with the result in reference ||. Note that the string susceptibility, ^stv, can 
be defined by (1) ~ e 2_ 7str > Thus we see that j s t r = — t for the k-th. multi-critical model. 

The one-point correlation function can be obtained by differentiating ( |38| ) with respect 
to t^. This gives for a general potential V(z), 

c 
Specializing to the k-th. multi-critical model gives 

r 1/k / / t \ e i+Ci+i)A 

M—l *">•{*-) = lk + m + ill + l)/k y («) 

Similarly the two-point function is 

(a h a h ) =-l— z h+l * log ( Z -^-\ = f U dx x h+l * = uh+h+1 (44) 

V h h ' v J 2vri S V z ) Jo (h+h + l) 

c 

Hence for the multi-critical model, 

e (h+te+l)/k 

<<^> = (, 1 + h + i) - < 45 > 

For the general potential we also have that u = (PP) V , a relationship which will prove 
useful later. The three-point function is 

/ dz z h+l2+h , A . 

(^ l2 a h ) v = f-^-jZy- (46) 

c 



d_ 

de 



c 



10 



(47) 



and in general, 



(cr h ai 2 ■■■ai s 



de 



dz j 

z L log 

2m & 



z-V 



C 



de 



u 



L+l 



(L + l) 



, (48) 



where L = J2 h- For the fc-th multi-critical model, we put u = e l ' k in this formula. 

»=l 
This result agrees with the calculation in references [g, || for the multi-critical models in 

topological gravity and also the similar result for the correlation functions in the multi- 
critical Hermitian matrix model [pH] ; this latter model is studied in section 0. 



3.3 Loop Operators 

In this section we define an operator W(l), which we will call a loop operator, by analogy 
with similar observables defined in matrix models. The form of (|42| ) suggests that we 
should define the loop operator as 



W ® = £ ~^n- 



n=0 



??.! 



(49) 



Let us now calculate expectation values of products of W(l) for a general potential and 
for the multi-critical models. The results in the previous section immediately give for the 
potential V(z), 



<"'(>i)>, ■#^ e,1 *(*- 1 lo 8( £ 7^ 



(50) 



c 
For the fc-th multi-critical model this gives 

(W(h)) = - / 
Jo 

For s > 2, 



dx e llX (x k — e] 



d xfc 



Jit 1 '" 



h 



(51) 



i[w(k 



\i=l 



d^ 


8-2 


[-/ 


dz 


Oe, 


) 


2vri 






L C 




d^ 


8-2 

) 


~e Lu - 


-l" 


de, 


L 





S explz^llo; 



j=l 



(52) 
(53) 
where L = ^ li as before. The fc-th multi-critical model is given by setting u = e ' in 

i=l 

this last equation. The relationship between these loop operators and the corresponding 



correlators which occur in matrix models is discussed in section 4.4. 
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3.4 The Dilaton and Puncture Equations 

The dilaton equation can be derived from (|38| ) in a simple fashion; we have 







Efc/sra 



k=0 



dti 



dz 
2vri 



V(z - V) log 



z-V 



(54) 



c 



and this gives us the dilaton equation, 






(1), 



0. 



(55) 



We will also derive the puncture equation. We have 

d 






(1>„ 



dz 
2^i 



c 



dz 



(z-V) 



(z - V) log 



z-V 



(56) 



After integrating by parts the right-hand side can be evaluated to give the puncture 

equation, 

I"/} S2. F^^ 1 „ 

(l)v = o e • ( 57 ) 



--^(fc + l)t fc+1 — 

k=o dtk 



de 



Many of the results in these sections have been derived by different methods in a 
number of papers [||, ||, [l(| . So far the contour integral representation of 2d topological 
gravity, at genus zero, has allowed us to reproduce the formulae for various correlation 
functions in a very simple fashion. However, the real power of the method will become 
apparent when we consider the Hermitian one-matrix model in the next section. 

4 Hermitian Matrix Model 

4.1 Contour Integral Representation 

Consider an Hermitian matrix model with an even potential; this is sometimes referred 
to as the reduced Hermitian matrix model. The partition function is defined as 



'■>N( 



DM exp 



-Tr U 



where M is an iV x N Hermitian matrix and the potential U is given by 



1 



U(M) = -M 2 + J2 



9p 



p=2 



NP- 



■M 2p . 



(58) 



(59) 



Note that we are writing the variables with bars on them in order to differentiate them 
from the rescaled variables, which will be introduced in the next section. The free energy 
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at genus zero is given by 



E (g) = lim 



1 , ( Z N {g) 



(60) 



The problem of calculating Eo(g) has been solved using the technique of orthogonal 



polynomials and in reference [21| it is shown that 



E (g) 



dr w Or) [1 — w(r)] lo 
o ' ~ V r 



w(r) 



where 



^ (2p)! p 



and a 2 is the solution of 
Integrating by parts gives 

Ml 



w{r 



1 = w(af 



p=2 



p!(p — 1)! 



a 2 >0. 



1 f a n ^2 d 
- dr (1 - lor-r- 

2 7o dr 



log 



w(r) 



(61) 
(62) 
(63) 

(64) 



We would like to rewrite this equation as a contour integral in such a way as to make the 
connection with 2d topological gravity clear. First let us change variables from the set of 
parameters {g^} to the set {t p } by defining, 



5p(2p)! 



v p_1 



for p > 2. 



(65) 



p\(p- 1)1 (i-h)p 

Note that the set of parameters {t p } has two more variables than the original set {~g~p~}, 
namely e and t\; this gives us some freedom in the choice of these last two variables. 
When it is convenient we will write e as to. In terms of these new parameters and using 
the variable z in place of r, w is given by 



1 



w(z) -l = -[ z '- V(z')} , 



where 



cz 



l-h 



and V(z) is defined as for topological gravity, 



p=0 



(66) 

(67) 

(68) 



The integral in (p4j) is along the real axis from zero to a 2 , where a 2 is given by the 
condition w{a 2 ) = 1. This corresponds to the range z' = to z' = u, where u is defined 
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by u = V(u). This is precisely where the cut occurs on the real axis for topological gravity 
with the potential V and suggests that we should rewrite (B3) as 



E (g) 



£[1-«M1' 



— \ ( w ^ 

dz V z 



log 



w{z) — 1 



(69) 



The contour encircles the cut from z = to z = a 2 , which has been introduced by adding 
the final logarithm in the integrand. Having rewritten Eq as a contour integral all that 
remains is to eliminate w(z) in favour of V(z), using 



— 1 f dz , 

E ° = -2?t*H {Z - Vr 

c 



_d b ( z-V + e 
dz \ z 



log 



z-V 



(70) 



where the contour C is around the cut [0, u] on the real axis. This gives the free energy at 
genus zero for the reduced Hermitian matrix model, and clearly the result is very similar 
to the genus zero partition function of topological gravity, equation (|3q). Equation ( |65|) 
relates the parameters in the two models. 



4.2 Rescaled Matrix Model 



The subsequent analysis will be simpler if we rescale the matrix in the matrix model, 
removing all the factors of N 1 ^ 1 ' and introducing a parameter in front of the quadratic 
term in the potential. We introduce a matrix <J> = M(l — £i)~2 N~? , set g p = g^(l — t\) p 
for p > 2 and define g\ = —\t\. Then after substituting into equations ( |58| ) and ( |59"| ) we 
have 

Z N {g) = fv$ exp[-iVTr [/($)] 

where 

1 

_< 

2 



P =i 



(71) 
(72) 



and the genus zero free energy is 



*M = Ji5.-Wf$> 



The contour integral representation is now, 
^ 1 , s 1 f dz , 

C 



d fz-V + e 
dz V z 



log 



z-V 



(73) 



(74) 



with V(z) defined by ( |6q ) as before. The relationship between the two sets of parameters 

is 



9 P (2p)l 
p\(p-l) 



-t p € 



p-1 



for p > 1. 



(75) 
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For a given C/($) we may freely choose e (with e / 0) and then all the other t p are 
determined. 

The relationship between the reduced Hermitian matrix model and topological gravity 
is simpler if we differentiate Eq with respect to £/. For I > 1, 



dE Q -Si,! 1 I dz 



dti 2(1 -ti) 2e 2 J 2ni 
c 



2z l - 1 (z-V)log(- 



z-V\ z l ~\z-V) 



z-V + e 



(76) 



which we gain after differentiating under the integral sign, integrating by parts and drop- 
ping analytic terms in the integrand. However the last term in the integrand cancels with 
the first term on the right-hand side, giving 

dE 1 



and hence for I > 1, 



3fa-l>v ( 77 ) 



^ «(i~ 1)1 7^*21 



<cr?_i) = lim — - '\ ,./• (TV $ M ). (78) 

This equation relates, at genus zero, the one-point correlation functions in topological 
gravity with correlation functions in the corresponding reduced Hermitian matrix model. 
Setting I = 1 gives 

< p >v=if 00 ^T( Tc * 8 >' (79) 

which shows that the correlation function for genus zero surfaces with a single puncture 
in topological gravity, is proportional to the partition function for triangulated surfaces 
with a marked link, with a suitable identification of the coupling constants in the two 
theories, that is, equation (|75|). 

In a similar fashion, one can show by differentiating Eq with respect to e that, 

1=1 

Equation (|7q) can be generalized to multi-point correlators. Let us define a set of 
operators, {£;}, in the matrix model as follows, 

El = Nel ~ ll ^m^ Tr (^ 21 ) for '- L (81) 



Then we have, by repeatedly differentiating (|78j) with respect to ti x , 

U-if[<ri) = Km (^) 2 (f[H k ) for/,>l, (82) 



1 = 2 / y U=l 



15 



where (• • -) conn denotes the connected part of the expectation value in the matrix model; 
that is, 

(S^eonn. = Pli^h) ~ Ph) Ph) (83) 



'conn 

and so on. 



Thus it is clear that the operator <j\ in topological gravity corresponds to the operator 
£/ in the matrix model, although the relationship is slightly different for the first operator. 
Note that this identification differs somewhat from the usual one in the literature ||, ^, 



22, 22], as many papers confuse 07 with an operator Oi, which is related to the l-th 



multi-critical model. This point will be discussed further in section M. 

4.3 Power Series Expansion 

In this section the connection between topological gravity and the Hermitian matrix model 
is used to write the genus zero free energy of the matrix model as a power series. It is 



simplest to begin with the version of the matrix model defined in section ^lj. Using the 
relationship, 

£b = £ -^log(l-ii), (84) 

and (|7^) we have 

Rewriting the derivative on the left-hand side in terms of {g^\ gives 

^ _ mm (l-tx) , 1 

i^P9p^^ = — 7^\ p ) v + o- ( 86 ) 



p=2 



dg P 



However, (P) _ can be written as a power series; using equation ([H]), after summing over 
no and differentiating with respect to to, we have 



2+J2(k-l)n k 

ni=o n 2 =o \j=i [l r ) 2+ E(Jb-l)n fc ! \*=i / 



tj \ to fc = 2 



k=2 



But, 



and so, 



00 ,ni / 00 \ / 00 \ -1-Vfcn fe 



n 1= ni " \fc=l / \fc=2 



(l"*l) 



00 00 



<*%=e E-n 



c 2 \ 'v 

?12=0 ri3=0 j=2 



t^ 1 



nj\ ui-hy 
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00 

fc=2 



2+ E(fe-iK)! 

fc=2 



Thus we have 

dE~n 



P9p W 



E 



CO 

k=2 



{thY 2 + £(fc-l)n fc ! p=2 

k=2 



n 



3£(2p)! \ n 



n p ! V p!(p- 1) 



(90) 



where we are summing each m (for i = 2,3,4, •• •) from to oo, but are excluding the 
case for which all the rii are zero. Integrating this equation yields 



E (g) = - E 



Y, kn k ) - 1 



k=2 



{mY (2+ J2(k-l)n k )\ p=2 

k=2 



n 



n P \ \ p!(p- 1)! 



(91) 



Note that for the special case in which only g 2 ^ 0, this reduces to 

^) = -Erwr (-12 5)*, 

^— n + 2 )!n! 

n=l v ' 

which is the well-known formula counting planar <3? 4 graphs [25]. Using (|84[) gives 



(92) 



Ma) - - e r k =L J \ n 

Kl' 2+ £(fc-l)n fc ! p =1 

\ k=2 J 



J_ / g p (2p)\ y 
n P \ \ p!(p- 1)! 



(93) 



where the sum is now over all values of ni,ri2, ■ ■ ■ from to oo, with them not all zero. 
This formula completely solves the problem of counting planar graphs, whose vertices all 
have even coordination numbers. 

4.4 Matrix Model Loop Correlators 

Let us define loop correlators for the reduced Hermitian matrix model (In]) as follows, 

i °° (ir <& 2fcl 



W(Pl) = T7 E 



and for s > 2, 



N ^ 2fci+l 

fc 1= o V\ 



oo (Tr $ 2fcl • • -Tr <3? 2fcs 

w(pi,---, p .) = jv'- 2 Yu — " 



(94) 



. 2*1+1 2fc s +l 

fc lr ..,fc s = l Pi P 5 



(95) 



In fact, we will only be interested in the genus zero part of these functions, which we will 
denote Wo(pi, • • • ,p s )', these are gained by taking the limit N — ► oo. 
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Using equation (|78|) we can write Wq (p\ ) in terms of expectation values for the corre- 
sponding topological model, 

1 ^ (2k x )\ (v kl -i) v 



W ^ Pl) pi + fc E k^h - 1)! j,2fe 1+ l e fc 1+ l • 



(96) 



This gives using (|42|) , 



c 
where we have defined z\ = jepf. Integrating by parts gives 



(97) 



W (pi 



1 
pie 



c 



2m \ z\) V z 



(98) 



where V'(z) = 4~. For example, in the case of the Gaussian model, U($>) = \& 2 , we 
have that V = e and integrating gives 



Wb(pi) = ( — 

\Pie/ Jo 



,, (l -f) "-* 



Pi 



(99) 



in agreement with the known result [24]. 
The general loop correlator is given by 



W (pi,---,p. 



d 



d 



where 



E 



dU( Ps ) dU( P2 ) 
1 i) ^ l (2Jfe)| 



W (pi 



E 



i a 



dU(p) £r[P 2k+1 dg k ^P 2k+1 kl(k - 1)! e k ~ l dt k 

Acting with this operator on the contour integral ([)?]) gives for s > 2, 



Wb(pi,---,p s ) 



where 



^(s** 



fife 
2~TTi 



c 



n i 



.1=1 



3 
2 \ 2 



de 



(100) 
(101) 



log 



z-U 



1 2 



(102) 



(103) 



Note that Woipi, ■ • • ,p s ) does not depend on the value of e; this parameter in the mapping 
between the matrix model and topological gravity can be freely chosen. It is convenient 
to choose e = 4. Then for s = 2 the contour integral can be converted to a real integral 
and calculated giving 



W (pi,p 2 ) 



1 1 



2{p\-plY 



2 \p\ — U 2 P 2 — u 
Pl\ -9 ^P2\\— -2£>lP2 



Pi — u 



p\-u 



(104) 
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where u is the solution of u = V(u) for e = 4 and we have assumed that pi,P2 > 0. 
For s > 3, 



w^,.., r .) = ^(f[' )f± 



\i=l ^1 



C 



Li=l 



n i-f 



Z \ 2 

2, 



<9e 



s-3 



1 



z-F 



, (105) 



which is being integrated around the pole at z = u. This can be easily evaluated to give 



1 d 



wb(Pi,-,P.)-(xFs:) Iw^n — 



M9u 



?/ 



(106) 



where M(u) = j (1 — V'(ii)) and again e = 4. These formulae, which were calculated for 
the reduced Hermitian matrix model, are identical in form to those for the loop correlators 
of the complex matrix model J2|, showing the close connection between these two models 
at genus zero. 

Next we take inverse Laplace transforms of Woipi, ■ ■ ■ ,p s ) with respect to each of the 
variables pf to get a new set of loop operators, 



i ~ C~ 5 



w o(h) = Jf E FTT-^T ( Tr *** 



f£> r (^ + I) 



(107) 



and for s > 2, 

oo 



n 



ni,—,n s 



=1 r(m + ±) r(n s + ± 



'Tr $ 2ni •••Tr $ 2 ™ s 



Taking inverse Laplace transforms of ( |102| ) gives 

r/2 



s-4 



Hence for s > 2, 



2vri 



n 



Vke zk 



,=i r(§; 



s-1 



s-2 



log 



z-y 



(108) 



(109) 



e=4 



r(|; 



9 



s-1 



V^M^) (W(h)---W(l s )) e=A , 



(110) 



where L = J2 h as usual and the expectation value of loop operators on the right-hand 

i=l 

side is evaluated for the potential V{z) with e = 4. The one loop correlator is 



w (h 



r(i: 



L J + JW'i)U 



fill) 
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These equations relate the matrix model loop correlation functions to the topological 



gravity loop operators which were studied in section |3.3| . Using the results in that section 
one can rewrite (|H0[) as 



iS-2 



w {h,- --,1s) 



r(|) 



d 



Vh 777 ^ Ur U- 



dL 



s-l 







Ik 



8-2 



e Lu -l 



L 



(112) 



£=4 



It is worth noting that in these formulae u is proportional to e. For example, consider the 



case of pure gravity, U(Q>) = ^ 2 + 92^ A - If we calculate the corresponding V(z) then 
the equation u = V(u) gives 



u 



24 5 2 



y/1 + 48g 2 - 1 



(113) 



This can be inserted into the previous formula to give the multi-loop correlator for pure 
gravity in the reduced Hermitian matrix model. 

5 Mult i- Critical Hermitian Matrix Models 
5.1 Definition of the Matrix Model 



In this section we consider the multi-critical Hermitian matrix models and compare them 

with the corresponding models in topological gravity. The matrix model will be defined 

as 

r r iv i 

(114) 



Z N = I V® exp 



iV 
— Tr U($) 



where the potential is 



1 



u($) = -$ 2 -J2T^u n ($) 



n=l 



and the potential U n ($) corresponds to that for the n-th multi-critical model, 

1 



£M$) 



-& 



1 



$ 2 --d» 4 , 

12 ' 






n\{p — 1)! 
(n-p)!(2p)! ' 



(115) 

(116) 
(117) 
(118) 
(119) 



these potentials were derived in reference pQ ]. 

The relationship between this matrix model and topological gravity is basically the 
same as in section fOL except that we have rescaled $ by a factor of e~2; that is, if for 
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the given £/(<£), the variables {g p } are defined by equation (ff|), then the parameters in 
the topological model are given by 



9 P (2p)l 
pKp~ !) 



^ = -J^W fa ^l. (120) 



The potential for the topological model corresponding to the above matrix model is 

oo 

V( z )=e-Y / T^[a-z) n -l] (121) 

n=l 

and the k-th multi-critical model is gained by choosing e 7^ 0, Jul = 1 and jujfc = — lj that 
is, 

y(z)=e + z + (l-z) fc -l = Vfc, (122) 

which we have denoted Vk- Note that this is different from the potential usually chosen in 
topological gravity for the fe-th multi-critical model, namely V = e+z—z, this potential, 
which has t\ = 1, would give g\ = —^ and hence t/(<3?) = g^ . Lacking a $ 2 term in 
the exponential, this matrix model does not correspond to the usual picture of surfaces 
made by gluing together polygons. This is the source of much of the confusion that has 
arisen when comparing the multi-critical matrix and topological models. 

5.2 Correlation Functions in the Multi-Critical Topological Model 

We begin by calculating the correlation functions for the topological model with the 
potential V& given in equation (|122j ). Then in the next section we consider the connection 
between these correlation functions and those in the matrix model. 

For the /c-th multi-critical model, u is given by u = Vk(u) and hence, 

w = l-(l-e)s =l-(Ae)^, (123) 



where we have defined Ae = 1 — e. 

The partition function for the fe-th multi-critical model is given by substituting Vk(z) 
into (pq), 



i/k 



{l)v * = \!o dX ( Ae - {1 - X)k f = ~\ I € d V^-y k )\ (124) 

where we have changed variables to y = 1 — x. Thus dropping the regular terms, 

this scales in the same fashion as (|4l|), which was calculated for the A:-th multi-critical 
model using the simpler potential. 
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Next we define a set of operators, {0{\, 

0< = -£(-l) n ('Vn, (126) 

n=0 KH/ 

which are chosen so that they scale in the correct fashion. The operator 0\ gives a factor 



of —(1 — z) in the contour integral and hence using (42), 

(0, K = - J^'Zy ,<(A £ - J) ~ [l + k ; k m + 1) [M»«»»\ (127) 

where we have dropped the regular terms in the last expression. 
Using (|48|), the general case (for s > 2) is given by 

UOk) -(-^) / W (128) 



w=l 



Me 



s-2 r i 



(^ + 1) 



(Ac) 



(L+l)/fe 



(129) 



where L = J2 k- 

Thus, up to regular terms, the scaling with respect to Ae of correlation functions of 
0\ for the potential, Vk = e + z + (1 — z) h — 1, is the same as the scaling with respect 
to e of correlation functions of ai using the simpler potential, V = e + z — z ; these were 



calculated in section 3.2. Note that some papers confuse 07 with O/, however they are in 



fact related by equation (|126| ). 

5.3 Correlation Functions in the Multi-Critical Matrix Model 

In reference [pH], Gross and Migdal study correlation functions for the fc-th multi-critical 
Hermitian matrix model at the spherical level. The correlators that they calculate scale 
in exactly the same fashion as those we derived in the previous section for topological 
gravity. It will be instructive to examine the connection between the two calculations. 

In their paper, the free energy of the matrix model is — log Zjy, the partition function 

for this model can be written as in equation ( |114| ). After dropping a number of terms, 

which they claim do not contribute to the critical behaviour, the free energy is given by 

F, where 

d 2 F 

■oF = /(*) ( 13 °) 

and f(t, {Hn}) is the relevant solution of 

00 

t = f" ~ E Mn/ n - (131) 

n=0 
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We will take the last two equations as being the definition of F and then later consider the 

connection between F and the free energy of the matrix model. The correlation functions 

are defined as 

;) ;) 

(132) 



9 d „, , ., 

-F(t,{»n}) 



{Mn=0} 



It should be noted that they introduced various powers of N in their definitions of the 
scaling variables: t, f(t) and so on. However these factors of N do not affect the following 
arguments and can safely be ignored. 

To make the connection with topological gravity we consider the potential, 

oo 

V(z) = - t + z + (l-z) k -Y / Mn(l - z) n . (133) 

n=0 

This is just equation ( |!2lD with the change of variables, 

W = fJ-n + S n ,i - 5 n ,k (134) 

and 

oo 

e = l-t-J2»n. (135) 

n=0 

The potential gives the fe-th multi-critical model perturbed by terms corresponding to the 
operators O n . Then the equation u = V{u) can be written, 

oo 

J2»n(l-u) n . (136) 

n=0 

{PP) in this topological model as was shown in 





t=(l-u) k 


Thus we see that / = 1 


— u, where u = ( 


section [3.2|. However ^ 


= — Mj:, hence, 




#F = f = 1 - 



» — < pp >= 1 -^< 1 >- < 137 » 

Integrating, we have that F = — (1), up to some regular terms which we will drop. Thus 
F, which is supposed to scale in the same way as the free energy of the matrix model, is 



in fact proportional to the partition function of the topological model defined in ( 133 ), 
up to some regular terms. Noting that t = Ae, when \x n = for all n, then, 

d d 



dnh djii 



-F 



= (O h ---O ls ) , (138) 

{/J.n=0} 



since deriving with respect to \x\ just pulls down an 0\ operator. Note also that the 
correlation function on the right-hand side is calculated in topological gravity for the /c-th 
multi-critical model (we calculated this in the previous section), it is not a matrix model 
correlator. Thus Gross and Migdal have been calculating topological correlation functions 
in their paper. The crucial question to answer is the following one: "If F in equation 
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( |132[) is replaced by the free energy of the matrix model, do the correlation functions still 
scale in the same fashion?" The answer is "Yes, provided that F is replaced by 6 2 Eq." 
The factor of e 2 in front of the free energy is necessary, as without it there are extra terms 
in the correlation functions, which would change the scaling behaviour if any of the li 
were greater than k. 

The free energy of the matrix model, Eq, is defined as in equation (|73|). We wish to 
calculate the correlation functions, which are derivatives with respect to fj,^ of —€ 2 Eq, 
evaluated at the A;-th multi-critical model (that is, [x n = for all n). However we know 



l\ _d_ 
n) dtr. 



(139) 



from (|126D that 

"' n=0 

Using equations ( |77| ) and (|3C|) one can calculate the effect of this operator on € 2 Eq; note 
that these equations were left unchanged under the rescaling of <J> by e~2 , which we 
performed earlier. After some algebra, we have 



_d_ 



- («•*) 



W=0} 



-elog k 



(Ae) 1 /* 

dy (Ae - y k ) 



y 



ky 



fc-i 



l-y l-y k 



The leading singular behaviour is given by 



_d_ 



(e 2 E ) 



■k 



-(Ae 



i+(/+i)/fc 



{/in=0} (i+k+m + i) 

as expected. 

The scaling of Eq can also be calculated from ( |141[) since, 



(Oi 



(140) 



(141) 



*(*» 



and hence as expected, 



e 2 E 



{)J.n=0} 



d 
dfi 



( £ ^o) 



{Mn=0 } (k + 1) 



(Ae) 1+ i 



k 1 



{m„=o} (2A: + l)(fc + 1 



-(Ae) 2+ i ~ 



-(1), 



(142) 



(143) 



In principle one could extend this calculation to multi-point correlation functions, 
however the algebra rapidly becomes tedious. In order to prove that derivatives of e 2 Eo 
have the same leading singular behaviour as derivatives of F, it is convenient to return to 



reference [21]. This paper gives an alternative form of equation (|6l] 

Eo~(g) = 



o 



1 J (1 M ( r o(x) 

ax (1 — x) log ' 



.t 



where ro(x) is defined by x = u>(ro). This leads to 



E o(g) = 7T lo S e 



1 

7 2 Jy, 



dx (e — x) \ogu(x) 



(144) 



(145) 
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where u(x) is the solution of V(u(x)) = u(x) for the potential (68) with to = x. Thus, 



d 2 /, n 3 



—2 (e 2 £ ) = - + log e - log u(e) 



(146) 



8e 2 V ~»J 2 
Dropping the first two terms, which are regular in Ae, we can differentiate with respect 

to (i h , 



de 2 



d 



(e 2 E ) 



(PPOh) 

(PP) 



(PPO h ) , 



(147) 



where the last expression uses the fact that, u(e) = (PP) = 1— (Ae) 1 ' and we are keeping 
only the leading non-analytic term. Integrating gives, up to regular and sub-leading terms, 

d 



9fJ>h 



e 2 E 



{Mn=0} 



This can easily be generalized to give 

d d 



dnh dm 



(e 2 E ) 



{Mn=0} 



(O h ) 



~ (O h ■ ■ ■ o, 



(148) 



(149) 



Thus as claimed earlier, the correlation functions for the reduced Hermitian matrix model 
scale in the same fashion as the topological correlation functions, for the fc-th multi-critical 
model. 

It should perhaps be noted at this point that if one rescales all the variables by various 
powers of iV as in reference |20|, then in the spherical limit TV — ► oo all the sub-leading 
terms in (149) vanish (this rescaling works provided we have dropped the regular terms). 
However, there is a simpler way of removing the sub- leading terms. From equation ( f77j ) 
we have that 











(e 2 E ) = (O t ) v , 



dm dm+iJ v > v 

which is true for a generic potential V(z). Differentiating repeatedly and evaluating for 
V k gives 

d d \ d d 



h+iJ dm 2 dm s v ' 



{Mn=0} 



$m x dfjL, 

This equation is exact and there are no sub-leading terms 

a 

dpi 



(o h --o ls 



(151) 



Using (|139|) and (|120|), one can write -£- in terms of derivatives with respect to the 



variables e and g p , 



_d_ 
dm 



EC- 1 ) 



P =i 



, l\(p-iy. d d_ 
(l- P )\(2 P )\dip~d~e- 



(152) 



In principle one can use this to operate on e 2 Eg and hence rewrite the correlation functions 
(151) as matrix model expectation values. The summation in (152) brings down a term 
proportional to Tr [/;($), but the derivative with respect to e is more complicated and 
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gives a number of terms including one proportional to Tr U(&) in the expectation value. 
As a result the multi-point correlators are not particularly simple. Calculating the one- 
point correlation function is however straightforward, 



(Oi) 



In fact, defining 











d/j,i dfn +1 



e 2 E ) = Jim - (TV [E/j($) - U l+1 ($)}} 



JV^oo N 



U t = Tr [[/j($)-[/ w ($)], 



U, = Tr 



(153) 

(154) 
(155) 

dfii dm dfii + i 
we have for the first few multi-point correlation functions in the spherical limit N — ► oo, 

1 



Ui($) - -U($) 



and 



d d 



d 



(o h o h ) v -- 
(o h o h o l3 ) v 



8 d <^F 



dfi h dm 2 

d d d 



u h u l2 



N 



U h 



., , , K e'E =[- (U h U h U k 
dfjt, h dfi h dfi h V / V e I \ * «mu. 



N 



(157) 
(158) 



and 



(O^O^ 



3^'4/y 



N 



u h u h u h u h 



+-(U h [U l2 U h + U h U h + U h U k 



(159) 



The correlation functions in topological gravity and the Hermitian matrix model can 
be related more simply at the cost of reintroducing sub-leading terms. If we define, 



O^Oi-O^, 



then in the spherical limit we have 



(O h ) 



N 



U h 



(160) 



(161) 



and in general, 

\ i=2 V V y \i=l / conn. 

This completes our study of the connection between the correlation functions in the 
reduced Hermitian matrix model and those in topological gravity. 
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6 Conclusion 

We have argued that it is possible in 2d topological gravity to consider observables de- 
pending on geodesic distances, such as the volume of space-time and the length of the 
boundaries. We then constructed a detailed mapping between topological gravity and the 
reduced Hermitian matrix model, both for operators and coupling constants. In particu- 
lar, the operators in topological gravity were associated with very specific matrix model 
operators even before the continuum limit of the matrix model was taken. As a spin-off 
of this identification we could use results from topological gravity to provide the complete 
solution of the counting problem for planar graphs whose vertices have all even coordina- 
tion numbers. It is an interesting problem to consider whether these counting results can 
be extended to arbitrary genus. Further, the simple and exact relation between a set of 
planar graphs and the partition function of topological gravity hints that there may be a 
simpler interpretation of the intersection indices on the moduli spaces of Riemann surfaces 
with punctures. However, a substantiation of such a claim requires first a generalization 
of the results of this paper to surfaces of arbitrary genus. An alternative representation 
of the operators a n might allow us to understand one remaining puzzle in topological 
gravity, namely the transition from local operators a n , which create only punctures in 
the Riemann surfaces, to operators which create macroscopic boundaries. Using matrix 
models such a transition is conceptionally simple since Tr $ 2n will create a macroscopic 
boundary, if n — > oo as the square root of the (discretized) volume. Since we have shown 
that (<r n _i) is closely related to Tr <£ 2n , the possibility of a geometric interpretation of 
(cr n _i) for n — ► oo should exist. We hope to be able to return to this question in a later 
publication. 
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